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Abstract 

We consider QCD instanton-induced contributions to lepton pair production in hadron- 
hadron collisions. We relate these contributions to those known from deep inelastic scattering 
and demonstrate that they can be calculated reliably for sufficiently large momentum trans- 
fer. We observe that the instanton contribution to the angular distribution of the lepton 
pairs at finite momentum transfer strongly violates the Lam- Tung relation - a relation be- 
tween coefficient functions of the angular distribution which is valid within the framework of 
ordinary perturbation theory. The drastic violation of this relation, as seen in experimental 
data, might be related to such instanton-induced effects. 
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1 Introduction 



The Standard Model of electroweak (Quantum Flavor Dynamics (QFD)) and strong (QCD) 
interactions is extraordinarily successful. This success is largely based on the possibility to apply 
ordinary perturbation theory to the calculation of hard, short- distance dominated scattering 
processes, since the relevant gauge couplings are small. Certain processes, however, can not be 
described by ordinary perturbation theory, no matter how small the gauge coupling is. These 
processes are associated with axial anomalies [1] and manifest themselves as anomalous violation 
of baryon plus lepton number (B + L) in QFD and chirality (Q 5 ) in QCD [2]. They are induced 
by topological fluctuations of the non-Abelian gauge fields, notably by instantons [3]. 

A number of non-perturbative issues in the Standard Model can be understood in terms of 
such topological fluctuations and the associated anomalous processes. On the one hand, QCD 
instantons seem to play an important role in various long-distance aspects of QCD, such as 
providing a possible solution to the axial U(l) problem [2] or being at work in chiral symmetry 
breaking [4]. In QFD, on the other hand, analogous topological fluctuations of the gauge fields 
and the associated B + L violating processes are very important at high temperatures [5] and 
have therefore a crucial impact on the evolution of the baryon and lepton asymmetries of the 
universe [6]. 

Are manifestations of such topological fluctuations also directly observable in high-energy scat- 
tering processes? This question has been seriously considered in the late 1980's, originally in 
the context of QFD [7]. But, despite considerable theoretical [8] and phenomenological [9] ef- 
forts, the actual size of the cross-sections in the relevant, tens of TeV energy regime was never 
established (for recent attempts, see Ref. [10]). Meanwhile, the focus switched to quite analo- 
gous QCD instanton-induced hard scattering processes in deep inelastic scattering [11], which 
are calculable from first principles within instanton-perturbation theory [12], yield sizeable rates 
for observable final state signatures in the fiducial regime of the latter [13,14], and are actively 
searched for at HERA [15]. Moreover, it has been argued that larger-size QCD instantons, be- 
yond the semiclassical, instanton-perturbative regime, may well be responsible for the bulk of 
inelastic hadronic processes and build up soft diffractive scattering [16]. It was emphasized for 
the first time in Ref. [17] that single photon or single W production at large transverse mo- 
mentum offers a possibility to study QCD instanton-induced effects from first principles at the 
LHC. Unlike the processes considered in the present paper, the dominant subprocess for this 
dedicated instanton search at high energies at the LHC [17, 18] is induced by gluon fusion, e.g. 
gg —>■ V + X, V = 7* —>■ £ + £~ . Moreover, the kinematical region is remarkable different from our 
region of interest, i.e. the available transverse momenta and virtualities are significantly larger 
than those we concentrate on throughout this paper. 

In this paper, we consider QCD instanton-induced contributions to lepton pair production in 
hadron-hadron collisions 1 (cf. Fig. Q). We relate these contributions to the ones previously 
calculated for deep inelastic scattering [12], thereby demonstrating that the former - like the 
latter - can be calculated from first principles. In particular, as already emphasized in Ref. [17], 
the typical inverse hard transverse momentum scale qj 1 in lepton pair production provides a 

1 This is often called the Drell-Yan process [19]. Instanton contributions to this process have been first discussed 
in Ref. [20] at a qualitative level. 
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Figure 1: QCD instanton-induced contribution to lepton pair production in hadron-hadron colli- 
sions, h\ + /t2 — ► (rif — 1) [c[r + qn] + £ + + t~ + n g g + X, corresponding to rif light flavours. 

dynamical infrared cutoff for the instanton size parameter p, thereby allowing for a controlled 
semiclassical approximation, which rests on the smallness of the QCD coupling at the effective 
momentum scale l/(p): a s {l/{p)) <C 1. Hence, in addition to deep inelastic scattering, lepton 
pair production at large transverse momentum may be viewed as a distinguished process for 
studying manifestations of QCD instantons. 

We put special emphasis on the angular distribution of the lepton pairs at finite momentum 
transfer. We observe that the instanton contribution strongly violates the Lam- Tung relation [21] 
between coefficient functions of the angular distribution, which has been verified within the frame- 
work of ordinary perturbation theory - the QCD improved parton model - up to 0{a 2 s ) [22,23] 
and even holds for the inclusion of parton transverse momentum and soft gluon effects [24,25]. 
Indeed, it has been argued that the drastic violation of this relation, as seen in experimental 
data [26-28], might be due to a non-trivial structure of the QCD vacuum [22], and in particular 
could be related to instanton-induced effects [29]. 

The outline of this paper is as follows: In Sect.|2]we introduce the instanton-induced contribution 
to lepton pair production. Afterwards we review the known results in the related process in DIS 
in an instanton background. The crucial part of this section is the continuation of these results to 
hadron collisions which leads us to the photon production tensor on partonic level. In Sect. El we 
will use these results to calculate the angular distribution of the produced leptons on the partonic 
(Sect. 13. 1|) and the hadronic (Sect. I3.2|) level. In Sect. 13.31 we give an outlook on the inclusion of 
multi gluon processes which lead to an enhancement of the instanton contributions. We present 
our conclusions in Sect. |U 

2 From deep inelastic scattering to lepton pair produc- 
tion 

We start with the derivation of the instanton-induced contribution to lepton pair production 
on the parton level. We will concentrate on the case with quarks in the initial state. These 
contributions dominate over the ones involving initial state gluons, at least for scattering processes 
where valence-like quarks and antiquarks contribute, e.g. in pp or -tr^N collisions. This is certainly 
different at very high energies where very small parton momentum fractions x dominate. Since 
the lower bound on x is set by M 2 /S, where M 2 is the invariant mass squared of the lepton 
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Figure 2: Instanton-induced process for n/ = 3, ul + ml — ► 7* + dn + dR + sr + sr + n g g, 
in leading semiclassical approximation. The amplitude involves the products of the appropriate 
classical fields (lines ending at blobs: fermionic zero modes (straight) and instanton gauge fields 
(curly)) as well as the non-zero mode quark propagator in the instanton background (quark line 
with central blob). 



pair and S is the hadron-hadron center of mass energy squared, the contributing values of x 
considered in our study are not so small for our chosen values of M 2 and S, see Section 13.21 For 
the main case of phenomenological interest, i.e. rif = 3 light flavours (m q (p) <C 1, for q = u,d,s), 
instanton-induced quark anti-quark annihilation involves in the final state at least two quarks 
and two anti-quarks of different flavour, such that the chirality is violated by 2n/ = 6 [2], plus an 
arbitrary number of gluons (g), e.g. (cf. Fig. Hj) 

u L + u L -> 7* + d R + d R + s R + s R + n g g . (1) 

t + r 



The amplitudes for the related processes in deep inelastic scattering, 

7* + g -»• u R + u R +d R + d R +s R + s R + (n g - 1) g, 



(2) 



have been derived, in leading-order semiclassical approximation, in Ref. [12]. For clarity and 
simplicity, let us concentrate here on the explicit result for the simplest appropriate 2 case n/ = 



n„ 



1 (cf. Fig. El (left)), 
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2 Note that the even simpler case n g = is not relevant for us, since the corresponding process in lepton pair pro- 
duction would contribute only at vanishingly small transverse momentum where, anyhow, instanton perturbation 
theory is not applicable (see below). 
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Figure 3: Instanton-induced process for rif = n g = 1, 7*(g) + g{p) — > qn(ki) + qR(k 2 ), in leading 
semiclassical approximation (left) and the analogous process, n f*(q)+g(p) — > q~R{ki) +^(^2), from 
ordinary perturbation theory (right). Both figures from Ref. [12]. 



and confront it with its chirality conserving counterpart from ordinary perturbation theory (cf. 
Fig. El (right)), 



(5) 



_ (9 - h) _ _ (<? 
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Here, e q is the quark charge in units of the electric charge e, g s is the strong coupling, A a , 
a = 1, ...,8, are the Gell-Mann SU(3) generators, and \i and p! are the four-vector indices of the 
photon and gluon, respectively (cf. Fig. EJ). The two-component Weyl-spinors Xl,r in Eqs. © 
and © satisfy the Weyl-equations, kxi(k) = 0, kxn(k) = 0, and the relations Xi{k) x\{k) = k, 
Xii(k) Xnik) = k. We used the abbreviations, k = k^a^, k = k^a* 1 , for any four- vector k^, with 
the familiar cr-matrices, cr^ = (1, a) and cf M = (1, — a), with a being the Pauli matrices. 

The other parameters in Eq. © arose from (the integration over) the instanton size distribution [2, 
30], whose two- loop renormalization group improved form [31], 
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(6) 



has been exploited, where p r is the renormalization scale and 

is the strong fine structure constant at two-loop, with 

A) = ll-~n /; ft = 102-^71/ (8) 
being the familiar perturbative coefficients of the QCD beta-function. The constant d is given by 
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Figure 4: Instanton-induced process for iif = n g = 1, qi{k\) + qi(k 2 ) — > 7*(?) + 9(p), in leading 
semiclassical approximation. 

with d = 0.466, C 2 = 1.51, and C 3 = 0.292, in the MS-scheme [32-34]. The variable b in Eq. © 
is a shorthand for the effective power of pp r in the instanton size distribution (|5jl. 

^/9o + ^^(A-12/5o). (10) 

47T 

It is important to note that the perturbative expression © for the size distribution is valid for 
small pA <C 1 where A is the fundamental scale in QCD. Indeed, a comparison with lattice data 
from quenched (n/ = 0) QCD [35] yields pA <0.4 for the fiducial region of instanton perturbation 
theory [36]. This can be translated into a fiducial kinematical region for instanton perturbation 
theory in deep inelastic scattering. On account of the fact that the main contribution to the 
integration over the instanton size comes from [12] 

6 + 3/2 b + 3/2 6 + 3/2 

corresponding to different terms in Eqs. (jHJ) and (|4*jl . one has to require that all virtualities, 
y/—q 2 , J—(q — k\) 2 , and J—(q — A; 2 ) 2 , exceed <2 m i n ~ (4 — 6) GeV, in order to stay in the realm 
of instanton perturbation theory. 

It is now straightforward to obtain the corresponding amplitudes relevant for lepton pair produc- 
tion via quark anti-quark annihilation, namely the one for the chirality violating instanton-induced 
process + <Zl(^2) - ► 7*(<?) + g(p) (cf. Fig. HJ) and the one for the analogous ordinary pertur- 

bative process (^(£4) + qR{k 2 ) — > 7*(<?) + g{p)- In fact, these processes are basically T-conjugates 
of the deep inelastic processes from Eqs. and (J3J), and the respective modulus-squared ampli- 
tudes are therefore identical, up to reflections of three-momenta. Some care has of course to be 
taken with respect to the photon virtuality: whereas in deep inelastic scattering it is space-like, 
Q 2 = —q 2 > 0, in lepton pair production it is time-like, M 2 = q 2 > 0. We will comment on this 
later. 

We have calculated the contribution of our simple processes to the partonic tensor 3 for inclusive 
quark anti-quark annihilation into a virtual photon, 

00 

w t , u {k 1 ,k 2 ;q) = ^w$(k u k2-,q), (12) 

n=0 
1 f 

"V ] (h, k 2 ; q) = — J dPS H T^ki^^q^x, ... ,p n )T*(ki,k 2 ,q,p u ... ,p n ), (13) 

3 Averaging over colour and spin of the initial state is implicitly understood in Eq. (|12(l ; the index n is to label 
besides the final state partons also their spin and colour degrees of freedom. 
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Following Ref. [21], this tensor can be decomposed as 



2 - q - pi 



Pr, 



(14) 
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— £ + kpkyWi, 



(15) 



with K' = ki + k 2 , k' = ki- k 2 , = g^ - q^v/q 2 , K'^ = g^K'^/y/s, and fc^ = g^k'" / y/s, 
where s = (&i + k 2 ) 2 . In turn, the different Wi can be obtained, 



i = l, 2,3,4, 



(16) 



from the partonic tensor with the help of the following projectors, 
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(17) 
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(19) 
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where t = (q — ki) 2 and u = (q — k 2 ) 2 . We are especially interested in the contribution w^J to 
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^-5^\p 2 ){2^Y^\k 1 + k 2 -q-p)w' l = U + \(k 1 + k 2 -q) 2 )w' i . (23) 
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Along these lines, we find that the simple instanton-induced process qi + q^ — > 7* + g contributes 
as follows to the functions w' ; , 
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We have obtained these results starting from Eq. contracting it with the gluon polarization 
vector e M (p) and taking the modulus squared, exploiting FORM [37] for the spinor traces. The re- 
sults (JH|)-(|2EJ) f° r the contribution of our simple instanton-induced process to the partonic tensor 
for inclusive quark anti-quark annihilation into a time-like photon look quite similar to the con- 
tribution of the analogous simple instanton-induced process to the deep inelastic structure tensor 
of a gluon found in Ref. [12]. A notable difference is the appearance of the factor Re(— l)^ 1 )/ 2 , 
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which reduces to unity in the space-like kinematics of deep inelastic scattering and was therefore 
not visible in the results of Ref. [12]. On the other hand, the full instanton contribution to the 
deep inelastic structure tensor of a gluon can be obtained from Eqs. (J2II)-(J2EI) by replacing the 



combinatorial factor l/((2 • N c ) 2 ) by 1/(2 • (iV c 2 - 1)) and by substituting Re(— 1 



,(6+l)/2 



by 1. 



As a check, let us also quote the corresponding perturbative contributions to the inclusive partonic 
tensor arising from Eq. (|5j). 
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We obtained this well known result, see e.g. [42], by exploiting the same FORM routines as the 
ones for the instanton-induced contribution, except for replacing the input amplitude (j3J) by the 
perturbative amplitude ©. 



3 The angular distribution of the lepton pairs 

In this Section, we will concentrate on the angular distribution of lepton pairs in instanton-induced 
processes - mainly concentrating on the simple one from the previous Section - and compare it 
to the one predicted from ordinary perturbation theory. 

In general, the angular distribution of the charged lepton £ + in lepton pair production, 

h 1 (K 1 ) + h 2 (K 2 ) -> l*(q)+X, (37) 

is described by three functions, A, /i, and v, which may depend on the kinematic variables of (|37|). 

' <h ° ' ( \ + A cos 2 6» + fi sin26> cos0 + - sin 2 6> cos20) , (38) 



a dQ 4tt A + 3 V 2 

9 and <j) being the polar and azimuthal angles of £ + , respectively [38]. These coefficient functions 
may be conveniently expressed in terms of hadronic helicity structure functions [39] , 

W T - W T 

A = TOT < 39 > 

" = WTwI' m 

V = 2W ^ . (41) 

Wt + W l k 1 
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Here, we exploit the so-called Collins-Soper frame [38], in which the frame dependent helicity 
structure functions read, in terms of the hadronic counterparts W{ of the previously introduced 
invariant functions to,-, 



where 
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(42) 
(43) 
(44) 
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(46) 



determines the transverse photon momentum q± with respect to the hadronic reaction plane. The 
kinematic variables S = (Ki + K 2 ) 2 , T = (q — K{) 2 and U = (q — K 2 ) 2 refer to the hadron level 
(cf. Eq. flHZD). Similarly, K = K x + K 2 and k = K x - K 2 . 



3.1 Parton level 

The contribution of our simple instanton-induced process as well as the contributions from or- 
dinary perturbation theory to these helicity structure functions are determined by folding their 
partonic counterparts with the parton density distributions. Before doing that it is instructive to 
consider first the partonic analogies of the quantities. 

Let us start with the contributions arising from ordinary perturbation theory, (cf. Eqs. (|3"T ]) -(|3l) |) ). 

7 ,,(P t )( 1 ) _ .MW O _ r 2 
\(pt)(l) = W T W L _ l ~ r U7) 

~ 4 pt)(1) + 4 pt)(1) 2 + 3r2' 140 

r/p t)(i) = ^i Pt)(1) _q ± 2s(s + M 2 )M 2 (t-u) 1 

~ ^ pt)(1) + w[ pt)(1) M(t-M 2 ) 2 + (u- M 2 ) 2 2M 2 s + ?,tu 1 ; 

2>(pt)(l) = - UQ) 

~ 4 pt)(1) + 4 pt){1) 2 + 3H- l4yj 

In terms of the partonic quantities, the ratio r 2 = q\/M 2 in Eq. ()46|) reduces to tu/(sM 2 ). In 
particular, we find the Lam- Tung relation [21], 

l_A(pt)(D _2z>( pt )« =0, (50) 

which, of course, holds also on the hadron level, 1 — — 2 = 0, as long as no intrinsic 
transverse momentum for the initial state quarks is invoked. It is nearly left intact even if one 
includes 0(a 2 ) corrections [23]. 
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The contributions arising from our simple instanton- induced process (cf. Eqs. (|24JI - (|28Jl ) are 
readily calculated along the same lines. They yield quite lengthy expressions, and we do not 
quote them all analytically, but will illustrate them, instead, graphically We stress, however, 
that the Lam- Tung relation is violated by instantons. This is apparent from the following, non- 
vanishing expression, 



" L ~ 2WAA ) ~ (t-M*){u-M*)[ { ' 




b+1 b+1 

sM 2 M 2 \ 2 M 2 \ 2 



tu \ —t J \ —u J J 
for the numerator of the Lam- Tung combination (cf. Eqs. ()39|) - (j41l) ). 



(51) 



l-\-2v= 2{WL - 2WAA) . (52) 

Note, that the factor (M 2 /\/tu) b+1 in the asymmetry (|5T)) arises from a non-planar diagram. 
That is in accordance with [29] where the importance of non-planar interference terms for the 
violation of the Lam- Tung relation were discussed. 

It is useful to view the partonic coefficient functions, for fixed M, as a function of r = q±/M and 
the partonic Feynman variable 

t — it 

x F = . (53) 

s 

In fact, A, fx, and v depend, for fixed M, only on r and x F . Their dependence on these kinematical 
variables is illustrated in Figs. El and El We observe the following features: 

(i) A (top panels) approaches —1 for large r for pure instanton- induced processes (dashed), i.e. 
these processes tend to be purely longitudinal (cf. Eq. ([39))) for large transverse momenta, in 
contrast to ordinary perturbative processes (dotted and Eq. (|4"7jl ). The total result for A (solid), 
taking into account both instanton and ordinary processes in the numerator and denominator 
of the partonic equivalent of Eq. ()39|) . shows little deviation from ordinary perturbation theory. 
Indeed, there are experimental hints for longitudinally polarized photons in hadron collisions 
towards larger xp [28]. Note that, even if the instanton- induced process qq — > j*g is suppressed, 
for larger xp, instanton effects might be relevant for this effect since gluon resummation leads to 
an enhancement in this kinematic region, see Section l3~B1 

(ii) The total result (solid) for fx (second panels from top) shows a quite significant deviation from 
ordinary perturbation theory (dotted) for sizeable xp and intermediate values of r. 

(iii) v (third panels from top) behaves quite differently in pure instanton-induced processes 
(dashed) and ordinary perturbative processes (dotted and Eq. ()49|l). In fact, instanton-induced 
processes have a value of 4 « 2 at small, but finite r and small xp , much larger as ordinary 
perturbative processes <^ 1). Also in the total result for v (solid) we observe a strong 

4 In fact, it follows from general arguments [22,29] that as long as instanton processes dominate over ordinary 
perturbative processes, one expects z>W w 2. 
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M = 7 GeV: 




g(0 10 r 

~~(pt) 5 



cr 







H i :: i I i 



I i N. 



I , + i I , In,, 1,1,^,1,1,1, S.^ 



Hhh-H- 



N 



0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 1 



q ± /M 



Figure 5: Coefficient functions of the angular distribution of lepton pairs from quark anti-quark 
annihilation, as a function of q±/M, for various values of xp, for M = 7 GeV. Dotted: result 
from pure ordinary perturbation theory; dashed: result from pure instanton perturbation theory; 
solid: total result from ordinary and instanton perturbation theory. Also shown, in the last row, 
is the ratio of the cross- sect ions of instanton perturbation theory to ordinary perturbation theory 
(dashed-dotted). In the numerical results shown we have chosen fj, r = M, A = 0.346 GeV, and 
n f = 3. 
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M=10 GeV: 
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Figure 6: Same as in Fig. but for M = 10 GeV. 
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enhancement at small, but finite r and small Xp in comparison to ordinary perturbation theory. 
Correspondingly, we find a strong violation of the Lam- Tung relation, which we display in the 
forth panels from top in terms of the parameter 

ft = _i(l-A-2i>) . (54) 

Whereas this parameter is identical zero in ordinary perturbative processes (dotted), it is about 
one, at small r, xp, for instanton-induced processes, leading to a drastic violation of the Lam- Tung 
relation in the total result (solid). 

(iv) Clearly, instanton effects in the coefficient functions are most visible in kinematical regions 
where the instanton induced cross-section a® dominates over the perturbative, a^ pt \ one. The 
instanton-induced features in Figs. El and El are indeed located where the ratio (bottom panels; 
dashed-dotted) 



r (pt) ~~ o„„(pt) , „„(pt) 



(55) 



becomes large. Obviously, it gets large towards small momentum transfer. The dominance of in- 
stantons is seen to set in like a "brick wall" . This sudden onset occurs practically at the boundary 
of the fiducial kinematical region of instanton perturbation theory, ~ Q m in or \J —u ~ Q m i n , 
for M > Q m in ~ (4 — 6) GeV. Therefore, the instanton features in the coefficient functions at very 
small momentum transfer, to the left of the sudden onset of instanton dominance in Figs. Eland El 
lie strictly speaking outside the range of validity of the semiclassical approximation. Fortunately, 
however, the coefficient functions are ratios of helicity structure functions (cf. Eqs. (j39j) - (}4*T]0 . 
and therefore the main uncertainties coming from the extrapolation of the perturbative expression 
of the instanton-size distribution cancel in them. Therefore, it is expected that our predictions of 
the coefficient functions remain also valid at smallish, but non-zero r. For very small q± = Mr, 
namely up to around 1 GeV, the simplest perturbative and instanton-induced sub-process qq — > 7* 
contributes and may change the angular distributions. 

(v) One feature of A^ pt ^ and z>( pt ) that translates unchanged even to the hadron level is the scaling 
behavior: they depend only on the ratio r. This is not the case when instanton effects are included 
since they vanish in comparison to the perturbative contributions for larger M 2 . This is basically 
triggered by the ratio cx® / which leads to a M 2 dependent weighting of the perturbative and 
instanton contribution. In addition, also A® and depend already slightly on M 2 . 

(vi) Remarkable is also the behavior for xf = 0. For vanishing q±, one expects to recover the 
well known leading-order angular distribution oc (1 + cos 2 ^), that is A = 1 and v, jl — 0. As 
one can see from Fig. El and El the function v tends towards 2. Note that v still vanishes in the 
limit q± — > for very small but finite xp- For large xp, the violation of the Lam- Tung relation is 
suppressed even in a region where the ratio a^ 1 ' / a^ pt > is not small. The very strong xp dependence 
will disappear after folding with the parton distributions as we will see in the next section. 

Strictly speaking, we should take n/ = 1 and a corresponding A = value for the calculation of 
the effective coupling parameters b, Eq. (fTU|). and Eq. (|2~5j). of the instanton contribution to the 
helicity structure functions, since our instanton-induced process corresponds to (the unrealistic 
case of) one massless flavour, the other flavours being integrated out. In the numerical results 
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Figure 7: Illustration of the dependence of our prediction of the Lam- Tung parameter k = 
t(1 — A — 20) on the choice of rif and A, for fixed M — 10 GeV and various values of xp. The 

thicker lines correspond to our default choice rif — 3 and A = A^ = 0.346 GeV, whereas the 

thinner ones correspond to rif — 1 and A = A^ = 0.241 GeV (other notations as in Figs. Eland 
EJ). 

shown in Figs. El and El we have chosen, instead, n/ = 3, and A = A^ = 0.346 GeV, for the 
calculation of the effective coupling parameters ^ pt \ and b. This value of A corresponds 
- according to the standard three- loop perturbative flavour reduction - to an n/ = 5 value 
^Ws = 0-219 GeV, leading to a running QCD coupling a s (m z ) = 0.119 at the Z-boson mass [14]. 
As illustrated in Fig. [7J our results for the coefficient functions are not largely affected if we 
choose instead the nominal value nt = 1 and a corresponding value for the A parameter, A^ = 
0.241 GeV. This value was obtained by a linear interpolation between the central values found in 
recent lattice investigations for n f = 0, A^ = 0.237 GeV [40], and n f = 2, A|| = 0.245 GeV [41]. 
Again, such details cancel to a great extend in the ratios of structure functions. This also refers 
to the dependence on the renormalization scale /z r , for which we have chosen M in the numerical 
results presented in Figs. El and El 
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3.2 Hadron level 



Since we have to deal with collisions of hadrons, the partonic Mandelstam variables s, t, u in Eqs. 
()24 )) -([28 )) are not observable. Firstly, we have to calculate the tensor (|12|) on the hadron level 
which involves a folding with the usual parton distributions, see e.g. Ref. [42], 

W»„(S,T,U) = y f--^E<(«.*.«) WsW+WM ( 56 ) 

3 J X\ X 2 i 

8n f dxi dx 2 r / s + t + u — M 2 \ „ . _ 
= y y — — ^ g J 2^w*(s,t,u) (qiMqifa) + qi(xi)qi(x 2 )) , (57) 

where the flavour dependence of is given by the relative charge e qi in £p~> (|2*9j) and ^ pt ^ flSTIj) . 
Note that the second equation (joTjl only holds for one parton in the final state, whereas the first 
equation is applicable for the general partonic tensor (|12|). The factors entering the hadronic 
tensor ()57|) are fixed in such a way that the tensor fits with the one defined in [42] 5 . We have to 
project the hadronic tensor (J57|) now on the hadron momenta K\ and K 2 to get the accessible 
hadronic structure functions Wi, 

W, u = -g^Wi + K,K V W 2 - K ^ + K ^ w 3 + kjUVt . (58) 

Here we have defined, similar to the partonic case discussed before, the vectors, K = K\ + K 2 , 
k = K\ — K 2 , = <7 M „ K v j\/S and = k u /y/S. Note the differences between the partonic 
momenta k', K' and the hadronic ones k, K in the hadronic tensor (|58|) . Due to the different 
projections on the hadron level the hadron structure function Wi is a linear combination of foldings 
of the four partonic functions w[ with the parton distributions. Using the partonic functions (J2*5j) 
- (j2Hj) for the instanton-induced contribution and Eqs. (J3*2~j) - for the perturbative one, we 
have now everything at hand to calculate the observable angular distributions (|3Uj) - (j41j) on the 
hadron level. 

For fixed M, the angular distributions on the partonic level depend only on r and xp, but they 
are independent of the center-of-mass (cm.) energy a/s, because the latter is fixed due to the 
relation s + 1 + u = M 2 . On the hadron level, however, the momentum fractions x\ and x 2 are 
variable, and the angular distributions depend, correspondingly, for fixed M, on the hadronic 
cm. energy y/S. Furthermore, the variable xp has to be replaced by the hadronic one, 

T — U 

Xf = — — , (59) 

which can be interpreted as the longitudinal photon-momentum fraction with respect to the 
momentum of the hadron h\. 

Figure |H1 shows the resulting angular distributions for proton-proton collisions at v^S* = 15 GeV 
for M = 7 GeV, Xp = 0, 0.1, 0.3, and varying values of r = q±/M. For the renormalization scale 

5 Constant factors are for our purposes actually not important since we are only interested in ratios of functions 
Wi. 
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we have chosen \x T = M and for the parton distributions the CTEQ6 dataset [43] 6 . 

The main difference to the partonic quantities, shown in Fig. 03 concerns the xp resp. Xp 
dependence. The strong xp dependence is smeared out on the hadron level. This smearing leads 
in particular to a suppression of the instanton- induced effect at small Xp. This has nothing to 
do with the parton distribution functions entering the hadronic angular distribution. Actually, 
since the angular distributions are ratios of two foldings with parton distribution functions, their 
dependence on the type of hadrons in the initial states is rather weak. The difference is just 
a consequence of the x\p dependence of the partonic Feynman variable xp that leads to the 
smeared out Xp behavior after integrating over 2^2 . For similar reasons, the M 2 dependence of 
the angular distribution in the instanton background is stronger than on the parton level, since 
for smaller ratios M 2 / s also smaller values of xp contribute, see Fig. |H1 and El but note that due 
to kinematical reasons only smaller values of r = q±/M are accessible for larger ratios M 2 /s. 

The fiducial region of instanton perturbation theory on the parton level yf—u, M > Q m in ~ 

(4 — 6) GeV can be mapped on the hadronic variables. One can check that theses relations are 
fulfilled for all x\ and xi for large enough values of M and q±, namely M > <2 m in and q± > Q m i n - 
Also on the hadron level our results should hold for even smaller q±, since the uncertainty towards 
smaller q± drops out in the ratios of the angular distributions, see the discussion in Sect. 13.11 



3.3 More partons in the final state 

As already mentioned in the introduction, the discussed instanton-induced sub-process qq — > 7*^ 
with only one gluon and no quarks in the final state is quite instructive since it contains already 
the basic non-trivial feature of the instanton-induced Drell-Yan process, namely the helicity flip 
of the quarks in the initial state, which is related to the chirality violation and is essentially 
responsible for the violation of the Lam- Tung relation [22,29]. But the rate of this asymmetry 
induced by instantons was certainly underestimated in the previous sections since it is well known 
that the resummation of the events with an arbitrary number of final-state gluons leads to a large 
enhancement which eats up at least partially the suppression of the instanton-induced process 
qq — > j*g. In addition, also the number of involved quarks in the subprocess is not realistic, see 
Sect. El 

A complete calculation of the angular distribution for this general instanton-induced process is 
beyond the scope of this paper and will be attempted in the future. Let us roughly sketch the 
general features of the complete process. Whereas in perturbative processes additional final- 
state gluons are certainly suppressed by an order of a s , every additional gluon in an instanton 
background leads to an enhancement of the order l/a s . Summing over all processes with an 
arbitrary number of gluons n g leads to an exponentiation of the inverse coupling constant. The 
resulting factor, combined with the tunneling factor, exp[— 47i/a s ] (cf. e.g. Eq. ()29|)). can be 
written as exp[— 4ii/a s F(x')]. Here, the Bjorken scaling variable x' = Q 2 /(Q 2 + M\) appears, 

6 Actually a consistent treatment of instanton induced effects requires also parton distributions including 
instanton-induced parton evolution. Since this modification enters the perturbative and instanton contribution in 
the same way this effect would change the angular distribution only in sub-leading order of instanton perturbation 
theory. 
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r = q ± /M 

Figure 8: The plot shows the angular structure functions similar to Fig. |S] but on hadron level. Therefore 
an integration over parton distributions, e.g. for the proton, is included. Due to the variable momentum 
fractions x\ and %2 one has to specify an additional kinetic variable, e.g. the cm. energy \[S. 
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gure 9: Same 



where Q is the relevant momentum transfer and Mx is the invariant mass of the produced partonic 
final state. The so-called holy-grail function F(x') [8] is normalized to one for x' = 1 and decreases 
towards smaller x' and therefore larger Mx ■ Let us mention that in the electroweak theory, where 
the coupling constant is much smaller, this mechanism is absolutely necessary for the process 
eventually becoming observable in the high energy limit [7]. 

For the process discussed in the present paper, Q 2 is given by the partonic quantities —t, —u or 
M 2 , whereas, in general, M x = {k\ + hi — q) 2 ■ Therefore, the integrands of the functions Wi 
(joT|) involve a factor exp[— 47r/a s F(x')}. For M\ = (x' = 1), and therefore also for the process 
with one final-state gluon, the factor exp[— 4n/a s } in Eq. ({2H} is recovered. For positive Xp, the 
smallest x' is given for Q 2 = —t, 

x , = -t = -t = -2M 2 /S + x 1 (^AM 2 /S(l + r 2 )+X 2 - X F ) 
' -t + M 2 x s + u-M 2 ~ 2x 2Xl - x 2 {^AM 2 /S{l + r 2 )+X 2 + X F ) 

It is easy to check that x' rises slightly with r and decreases towards the largest accessible values 
of Xp. Therefore, we can conclude that the instanton-induced effect in v and k (see Figs. |HJ 
IS] and El) will shift to slightly smaller r. Furthermore, we expect a significant enhancement of 
the instanton effect for larger Xp. Correspondingly, the suppression of the simplest instanton- 
induced process at large Xp, which we observed before, might be compensated. Note that the 
applicability of instanton perturbation-theory now requires in addition a cut x'>x cnt , where x cut 
is approximately 0.35, see Ref. [36]. One can check that this requirement can be fulfilled for all 
Xi, x 2 and r as long as Xp is not too large, or for all Xp for large enough ratios M 2 /S. 

Beside the discussed instanton-induced multi-gluon process also other perturbative processes may 
contribute. Firstly, we have not taken into account an enhancement of perturbative contributions 
due to soft gluon resummation [24, 25] at small q± since the instanton-induced contribution that 
we have calculated is not reliable in this region anyway. In higher order a s also new processes 
contribute to the angular distribution which lead to a small violation of the Lam- Tung relation 
already in the purely perturbative framework [23]. For small tranverse momenta q± the usual 
factorization is not reliable anymore and transverse parton momentum distributions become im- 
portant. However, as already mentioned in the introduction, higher order contributions, soft gluon 
effects and parton transverse momentum are not able to explain the observed strong violation of 
the Lam- Tung relation [23-25] . 



4 Conclusions 

We have calculated the angular distribution of the produced leptons in hadron hadron collision in 
an instanton background. It turns out that, for large enough photon virtualities M 2 and transverse 
photon momenta squared q\ , only small instantons contribute. Therefore, the instanton-induced 
contribution is fiducially calculable in this kinematic region using techniques of instanton pertur- 
bation theory. The most remarkable property of the resulting angular distribution is the violation 
of the Lam- Tung relation which is conserved to very high accuracy in usual perturbation the- 
ory, but violated in experiments. This effect is a direct consequence of chirality violation in the 
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background of QCD instantons which leads to a non-trivial spin-density of the quark-antiquark 
pair in the initial state as it has been argued in [29]. Therefore, lepton pair production in hadron 
collisions is potentially a very good testing ground for instanton- induced processes: the violation 
of the Lam- Tung relation is reliably calculable in instanton perturbation theory and absent in 
usual perturbation theory. 

We restricted ourselves to the simplest partonic subprocess qq 7* + g. Since the inclusion of the 
more realistic general processes qq —> 7* + (rif — l)qq + n g g was beyond the scope of this paper, 
we cannot compare our results directly with the available data. However, the small violation of 
the Lam- Tung relation on the hadron level arising from the simplest partonic process is already 
quite promising, notably in view of the expectation that additional gluons lead to a substantial 
enhancement of the instanton-induced effect, as known from analyses of the related processes in 
deep-inelastic scattering and from the general arguments presented in this paper. 

Finally, let us mention that, beside further theoretical efforts, more experimental data are re- 
quired for testing instantons in the angular distribution of produced leptons at a hadron collider. 
Fortunately, there are new medium energy projects under way that are also dedicated to study 
the Drell-Yan process, e.g. at the forthcoming facilities GSI-FAIR [44] and J-PARC [45]. 7 Exper- 
iments at RHIC may also give further information on lepton pair production. In general, it seems 
that fixed target experiments are especially well suited for our purposes, since on the one hand 
the involved momenta are smaller and on the other hand the luminosities are larger. Therefore, a 
huge amount of lepton pairs should be observable which is absolutely necessary for reconstructing 
a whole angular distribution. 
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